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ABSTRACT 

Context. Circumstellar envelopes (CSEs) of a variety of evolved stars have been found to contain ammonia (NH3) in amounts that exceed predic- 
tions from conventional chemical models by many orders of magnitude. 

Aims. The observations reported here were performed in order to better constrain the NH3 abundance in the CSEs of four, quite diverse, oxygen- 
rich stars using the NH3 ortho J^ = \q- Oq ground-state line. 

Methods. We used the Heterodyne Instrument for the Far Infrared aboard Herschel to observe the NH3 7^: = lo - Oo transition near 572.5 GHz, 
simultaneously with the ortho-HiO Jk,„Kc = 1 1,0 - lo.i transition, toward VY CMa, OH 26.5+0.6, 1RC+ 10420, and IK Tau. We conducted non-LTE 
radiative transfer modeling with the goal to derive the NH3 abundance in these objects' CSEs. For the last two stars. Very Large Array imaging of 
NH3 radio-wavelength inversion lines were used to provide further constraints, particularly on the spatial extent of the NH3-emitting regions. 
Results. We find remarkably strong NH3 emission in all of our objects with the NH3 line intensities rivaling those of the ground state H2O line. 
The NH3 abundances relative to H2 are very high and range from 2 x 10"' to 3 x 10"' for the objects we have studied. 

Conclusions. Our observations confirm and even deepen the circumstellar NH3 enigma. While our radiative transfer modeling does not yield 
satisfactory fits to the observed line profiles, it does lead to abundance estimates that confirm the very high values found in earlier studies. New 
ways to tackle this mystery will include further Herschel observations of more NH3 lines and imaging with the Expanded Very Large Array. 
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^. 1. Introduction 



Ammonia (NH3) was the first polyatom i c mol ecule detected 
in an astronomical object (ICheung et al.l Il968h . It is ubiqui- 
tous in dark, dense interstellar cloud cores and an eminently 
usefu l thermometer of t hese regions (IWalmslev & UngerechtsI 
I1983L lDanbvetal.1ll988 l). This, aided by the easy observability 
of its inversion lines - many of the astronomically most impor- 
tant ones crowd around 1.3 cm wavelength (24 GHz frequency) 
- make N H^ one of the most frequently observed interstellar 
molecules (iHo & Towneslll983h . 

NH3 has also been detected toward a still limited, but di- 
verse number of CSEs around evolved stars, first using infrared 
(IR) heterodyne absorption spectroscopy toward the high mass- 
loss asymptotic giant branch (AGB), extreme carbon star CW 
Leo (= IRC+10216; iBetz etal.iri979l) . In addition, absorption 
was found toward a number of oxygen-rich objects that included 
the long-period variable (LPV) o Ceti (.Betz & Goldhabei.1985,) 
and the supe r- or even hyper-luminous objects VY CM a and 
IRC+10420 (iMcLaren & Betjll980i iMonnier et al.ll2000h . The 
last study finds that around VY CMa the NH3 is forming at > 40 
stellar radii away from the star, where dust formation has well 
started. 

Contemporaneously, several r adio inversion line s were de- 
tected also toward IRC+10216 dKwok et all IT98T1 iBell et all 
1 1982, Nguven-0-Rieu et al. 198-^. Later on, high-velocity cm- 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



wavelength NH3 emission plus absorption was found toward 
the bipol ar protoplanetary nebulae (PPNe) CRL 2688 and 
CRL 618 (lTruong-Bache"tanil988llMartin-Pintado & Bachilleii 
ll992ilTruonp-Bach et alj|1996h . Toward CRL 618, P Cygni pro- 
files are observed with a full width at zero power (FWZP) of 
w 100 km s"'. Menten & Alcolea ( 1995) detected high-velocity 
NH3 radio emission toward IRC+10420, the high mass-loss rate 
LPV IK Tau, and the PPN OH 231.8+4.2. In the last case, they 
find high- veloci t y emis sion over a; 70 km s"' FWZP. Recently, 
Has egawa et al.l ( l2006l) report and discuss observations of the 
NH3 line central to the present study, the ortho-NH3 lo - Oo tran- 
sition, toward IRC+10216 made with the Odin satellite. 

One common, surprising result of all the above studies is 
the exceedingly high NH3 abundances they report. Most of them 
cite values of several times 10"' or even 10"^ relative to molec- 
ular hydrogen. These numbers are in stark contrast to the re- 
sults of thermodynamical equilibrium calculations for the atmo- 
spheres of cool stars, which predict the production of only neg- 
ligible amounts of NH3, of order 10"'^ (.Tsuil.1964.) . The pio- 
neering stu dy of Tsuji has been confirmed by more recent work 
(see, e.g., iLafont et al.lll982l ICherchneff & Barkeilll992 [ both 
for C-rich CSEs). Somewhat ad hoc approaches to bringing ob- 
servations and theory closer together involved injecting a sig- 
nificant a mount of NH3 in the inner envel ope dNei ad & Millar! 
1988, Ne rcessian et al.lll989l) . IWillacv & C herchneffl (Il998h ^ 
elude shock chemistry in their model of IRC+10216, but still 
only produce an abundance of 4 x 10"'', at least three orders of 
magnitude below the value implied by observations. As to bona 
fide shocked regions like PPN outflows. iMorris et al.l(ll987i) sug- 



K. M. Menten et al.: Ammonia in 0-rich evolved stars 



gested that, for OH 231.8+4.2, N2, which binds most of the ni- 
trogen, might be dissociated in the high-velocity gas and that the 
high NH3 abundance might be the result of a series of hydro- 
genation reactions. Whether this can be confirmed by detailed 
chemical models remains to be explored. 

For the present study, as described in Sect. |2] we observed 
the NH3 Jk - lo - Oo ortho ground-state transition in O-rich 
stars of widely diff'erent natures and mass loss rates: the high 
mass-loss LPV IK Tau, the peculiar red supergiant VY CMa, 
the archetypical OH/IR star OH 26.5+0.6, and the hypergiant 
IRC+ 10420. We chose a receiver setting that allowed simultane- 
ous observations of the Jk,„Kc - li.o - lo,i transition of ortho- 
H2O. All of these objects have dense CSEs, and NH3 has been 
previously detected toward all of them but OH 26.5+0.6. In par- 
ticular, for IK Tau and IRC+10420, single-dish observations of 
the (J,K) = (1, 1) and (2, 2 ) inversion lines have been reported 
by iMenten & Alcoleal (1995). Moreover, the emission in these 
lines has subsequently been imaged with the NRAO Very Large 
Array (VLA) with a resolution of a few arcseconds (Menten et 
al., in prep.; see Sect. 12.31 ). 

The critical density of the cm-wavelength inversion lines is 
on the order of lO'* cm"-', while that of the sub-mm Iq - Oo 
transition has a value ~ 4 orders of magnitude higher Thus, 
both types of lines should provide complementary information 
on different regions of the envelope. In Sect. 13.11 we give a gen- 
eral description of the sub-mm spectra we obtained with HIFI 
aboard He rschel (Pilbratt et al. 2010, in this volume). Thereafter, 
in Sect. 13.21 we present radiative transfer calculations conducted 
to model the observed line profiles, taking advantage of the con- 
straints from the VLA imaging. These lead to NH3 abundance 
determinations. 



2. Observations 

2. 1 . Herschel/HIFI submillimeter observations 

The observations were made with the two orthogonal HIFI re- 
ceivers available for each band, which in all cases work in dou- 
ble side-band (DSB) mode (see de Graauw et al. 2010, in this 
volume). This effectively doubles the instantaneous intermedi- 
ate frequency (IF) coverage. We observed the four stars de- 
scribed above with a tuning that, in the upper sideband, cov- 
ers the frequency of the J^ - lo - Oo ground state transition 
of ortho-NH3 at 572.4981 GHz. The tuning was chosen to also 
cover the frequency of the Jk„,k,. = 1 1,0 - 1 0,1 ground state 
line of ortho-HaO at 556.9360 GHz in HIFI's lower sideband. 
The observations were obtained using the dual-beam-switching 
(DBS) mode. In this mode, the HIFI internal steering mirror 
chops between the source position and a position believed to 
be free of emission, which was certainly the case for our ob- 
servations. The telescope then alternately locates the source in 
either of the chopped beams, providing a double-difference cal- 
ibration scheme, which allows a more efficient cancellation of 
the residual standing waves in the spect ra. Additional details 
on this observing mode can be found in iHelmich et al.l (l20I0i 
in this volume). The double sideband system temperature was 
^ 100 K, and the calibration uncertainty is estimated to be 10%. 
Spectral baselines were excellent. HerscheVs, beam had a size of 
37" FWHM at the observing frequency, which is much larger 
than the NH3 -emitting regions of all our sources. 

The HIFI data shown here were obtained using the wide- 
band spectrometer (WBS), which is an acousto-optical spec- 
trometer, providing a simultaneous coverage of the full instan- 
taneous IF band in the two available orthogonal receivers, with a 



(oversampled) channel spacing of 0.5 MHz (0.27 km s '), about 
half the effective resolution. Spectra in the figures have been re- 
sampled and smoothed to a channel spacing of a; 1.1 km s ' . 

The data were processed with the standard HIFI pipeline us- 
ing HIPE, and nonstitched Level-2 data were exported using the 
HiClass tool available in HIPE. Further processing, i.e. blanking 
spurious signals, first order polynomial baseline removal, stitch- 
ing of the spectrometer subbands and averaging, was performed 
in CLASS. Since the quality of the spectra measured in both hor- 
izontal and vertical polarization was good, these were averaged 
to lower the final noise in the spectrum. This approach is justified 
since polarization is not a concern for the presented molecular- 
line analysis. All HIFI data were originally calibrated in units of 
antenna temperature (rp and were converted to the main-beam 
temperature (Tmb) scale according to Tmb - T^/tjmb, with the 
main-beam efficiency tjmb - 0.68. In all cases we have assumed 
a side-band gain ratio of one. 

2.2. Ammonia spectroscopy and astropiiysics 

The main focus of this letter is on the Jk - lo - Oo line of ortho- 
NH3 ■ Ammonia microwave spect r oscopy ha s a long history (see, 
e.g., iTownes & Schawlow''1955 UKukolich ' 1967). Very briefly, 
because of the possible orientations of the hydrogen spins, two 
different species of NH3 exist that do not interconvert, ortho- 
NH3 and para-NH3. Ortho-NH3 assumes states, Jk, with K - Q 
or 3«, where n is an integer (all H spins parallel) , whereas K + 1> 
for para-NH3 (not all H spins parallel). The principal quantum 
numbers J and K correspond to the total angular momentum and 
its projection on the symmetry axis of the pyramidal molecule. 

The temperature corresponding to the energy of the lowest 
para level (Jk - 1 1 ) is 22 K above that of the lowest ortho level 
{Jk - Oo). Therefore, for formation in the interstellar gas phase, 
which involves reactions with high exothermicities, the ortho- 
to par a-NH3 ratio is expecte d to attain its equilibrium value of 
unity (Umemoto et al. 1999). This situation is also expected to 
hold for CSEs, given that the IR studies cited above place the 
NH3 they observe in the hot medium close to the star. 

A high-resolution study of the NH3 Ip - Oo transition has 
very recently been presented bv lCazzoIi et al.l (12009). (Only) its 
upper state is split into several hyperfine structure (hfs) compo- 
nents with the ~ 1-2 MHz splitting resulting from the coupling 
of the quadrupole moment of the N nucleus with the electric 
field of the electrons. Two of these components are further split 
by magnetic interactions. The mean frequency is 572498.1 MHz 
and the centroid frequencies of the three main hfs groups are all 
within 2 MHz, corresponding to x; 1 km s ', much less than the 
line widths observed for the targets of this study (see Fig.[T]i. 

2.3. VLA observations of inversion lines 



In Sect. 13.21 we use data of the (J,K) = (1, 1) and (2,2) inver- 
sion lines to constrain our models for IK Tau and IRC+10420. 
For these stars, single-dish observations of those lines made 
with the Effelsber g 100 meter telescope were reported by 
Menten & Alcolea ( Il995h . In addition, to place constraints on 
the spatial distribution of the NH3 molecules, we have used data 
obtained with the VLA that will be published separately (Menten 
et al. in prep.). The hfs splitting in the inversion lines in velocity 
units is much wider than for the rotation line. However, because 
both the (1,1) and (2,2) lines are very optically thin, as indicated 
by the spectra and supported by our modeling (see Sect. I3.2l i. 
any contribution of the hfs components will be factors of several 



K. M. Menten et al.: Ammonia in 0-rich evolved stars 



weaker than the main hfs component and neglected in the mod- 
eling. The intensities of the spectra produced from the VLA im- 
ages, which were restored with a circular beam of 3. "7 FWHM, 
used in that section, are consistent with the published 100 m tele- 
scope values. 

3. Results and analysis 

3.1. NHj, versus H2O emission 



In Figs. [T] lA.ll and Table [T] we present the results of our HIFI 
observations of both the NH3 and the H2O ortho ground state 
lines, together with the results of our NH3 modeling. All our ob- 
served positions agree to within 2" with the stars' 2MASS posi- 
tions, whichthemselves have an absolute accuracy of better than 
0."1 (ICutri et al.l l2b03). The determination of the LSR ranges is 
somewhat subjective and the upper and lower velocities are un- 
certain by ~ a few km s"' for weaker lines. For all entries, the 
formal error in J Tusdv is smaller than 0.1 K km s '. For VY 
CMa we used the higher of the literature mass-loss rate values 
scaled to the recently measured trigonometric parallax distance, 
1 100 pc (Ch oietalJl2008l) . For both IK Tau and IRC+10420 the 
lower and the higher values of Xnh, are implied by the cm lines 
and the submm line, respectively. 

Inspecting Fig. [1] and the table, it is striking to see that the 
luminosity (integrated intensity) in the NH3 and H2O ground- 
state lines is of comparable magnitude for all of our objects. The 
NH3/H2O line ratios are ^ 0.50, 0.30, 0.28, and 0.42 for IK Tau, 
VY CMa, OH 26.5+0.6, and IRC-H10420, respectively. One has 
to keep in mind that H2O is a major molecular constituent of our 
CSEs, while even the presence of observable NH3 emission is 
completely unexplained! 

Another remarkable result is that the velocity ranges covered 
by the two lines are almost identical, which suggests that the 
bulk of the material producing the emission for both is similar. 
Moreover, for all our targets, both lines' FWZP values are lower, 
but comparable to twice the terminal velocity, implying that both 
molecules are present in the outer layers of the envelope, where 
the material has almost been fully accelerated. Furthermore, we 
point out the clear self absorption in the blue wing of the H2O 
lines toward IK Tau and VYCMa, which proves that the line 
emitting region covers the whole envelope. Whether this is also 
true for the NH3 line is a priori not clear 
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Fig. 1. HIFI spectra compared to radiative transfer model results: 
The thick black lines show spectra of the Jk = lo - Oo transi- 
tion of ortho-NH3 for (top to bottom) IK Tau, VY CMa, OH 
26.5+0.6, and IRC+10420. The blue Hnes represent the ortho- 
H2O Jk„,k, - li.o - lo.i spectra for the same stars in the same 
order The dashed red lines represent the predictions for the NH3 
line resulting from our radiative transfer modeling. The left- and 
righthand ordinates give the main-beam brightness temperature 
scales for the NH3 and the H2O lines, respectively. Except for 
IRC+10420, they have different ranges. 



3.2. Radiative transfer modeling and constraints on 
abundances 

The NH3 emission of the sources has been modeled with 
the Monte Carlo radiative transfer code RATRAN developed 
and described by iHogerheiide & van der TakI (2000|). For NH3 , 
RATRAN uses collision rates calculated bv lDanbv et all (11988^ . 
Power laws for the density and temperature were used to de- 
scribe the physical structure of the envelope, using as input 
published values for the mass loss rate and expansion veloc- 
ity (see Table [B- For IRC +10420, envelope parameters from 
iDinh-V-Trung et alj ( l2009h were used and for IK Tau, VY CMa, 
and OH26.5+0.6 w e refer to the modeling o f De cin et al. (2010), 
iDecin et all (l2006l) . and iJusttanont et alJ (12006), respectively. 
The VLA data yield an extent of ^ 4" for the NH3 emitting 
region around IRC+10420 and x 2."6 for that around IK Tau, 
numbers we use for our modeling. 

To fit the submm NH3 lines observed with HIFI, the NH3 
abundance was varied in a first iteration. For the two sources 
with additional data from the cm inversion lines, those (para) 



lines were modeled as well, using an ortho-to-para ratio of 1, 
appropriate for formation of NH3 under high temperatures (see 
Sect. 12.2b . Interestingly, this does not lead to a satisfying fit for 
both the cm NH3 inversion lines and the submm ground state 
line. With a fit adjusted to reproduce the cm lines, the submm 
line is underestimated by a factor of 10 for IK Tau. 

Since high densities are needed to excite the submm line, its 
emission must arise from the inner part of the envelope, further 
inward than the cm-line emitting region. Our modeling suggests 
densities above a few times 10^ cm "* and temperatures in the 
10-100 K range. This discrepancy between the physical condi- 
tions required to produce the observed cm emission, on the one 
hand, and the submm emission, on the other, is even greater for 
the very extended but relatively low-density shell of IRC+10420. 
To reach densities high enough to excite the NH3 sub-mm line, 
the mass-loss rate had to be increased to 2 x 1 0'-^ Mp yr~', i.e., 
a facto r 3 higher than the value derived by Di nh-V-Trung et aH 
(|2009|) and the inner radius reduced to 5 x 1 ' ^ cm, which is twice 
the value of the hot inner shell proposed by these authors. Then 
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Table 1. Results of HIFI NH3 and H2O observations and NH3 modeling 



Object Q'j20oo Snooo D M^' v2 olsr range j Tundv W^ r^, 9^, Xnh, 

(kpc) (Moyr') (Iflns-') (kms"') (Kkms"') (M0 yr') (cm) (") 

IKTau 03''53"'28.'8 +ir24'23" 0.2511) 8 x 10-*{2j 18 {2| [11.5,57.2] 5.1 8 x 10^' 5.0 x IO'M.3 (1-3) 10^*^ 

[10.4, 59.5] 10.3 

VYCMa 07 22 58.3 -25 46 03 1.1(31 (0.8-3.4) 10-''{4) 35 |4j [-23.7,64.3] 15.5 1.8 x lO^* 3.5 x 10"^ 2.1 3 10"^ 

[-34,79] 51.2 

OH26.5+0.6 18 37 32.5 -05 23 59 1.37 |5| 2.4xl0-*(5) 16 16| [14.6,43.2] 0.96 2.4 x lO"'* 2.5 x 10'« 1.2 3 10^' 

[15.0,43.2] 3.4 
IRC+10420 192648.1 +112117 5 (7) 8x10-^7) 38 17| [39,112] 3.1 2 x lO-"* 1.5 x 10" 2.0 (2-5) 10"' 
[42.8,116.4] 7.4 

Notes. Columns are (from left to right) the object, J2000 right ascension and declination of observed position, distance, mass-loss rate from 
the literature, terminal expansion velocity (from CO data), FWZP LSR velocity range with observed emission, integrated main-beam brightness 
temperature, as well as mass-loss, outer radius of the NH3 distribution (in cm and in arcseconds), and [NH3/H2] abundance ratio derived from our 
modeling. For the dlsr range and J Tusdv for each source, the first and second rows give the values for th e NH3 and the H 2 O lin es, r espectively. 
Numbers in braces relate to the following references for distances and mass loss rates: {1) Hale et al. ( 1997), |2| Deci n et aLU201(]h . l31 IChoi et al.l 
(.2008,) . {4i iDecinetaL(,2006) . i5I Justtanont et a l.,(2006). (6l iKemper et al..l,2003) . !7l ,Dinh-V.-Trung et al., (2009). 



agreement between the abundances obtained from the cm and 
submm lines can be reached within a factor of 2. Interestingly, 
De Beck et al. (2010, accepted for publication), derive a value 
of 3.6 X lO""* Mq yr"' for IRC+10420's mass loss rate based on 
multi-transition modeling of CO. 

To reproduce the strong submm NH3 line from VY CMa, 
the highest mass-l oss rate and largest radius from the various 
shells discussed bv lDecin et al.l |2006) had to be used, scaled to 
D = 1 100 pc (see Sect. 13.11 ). The outer radius that led to a best 
fit for OH26.5-H0.6 is 2.5 x 10'^ cm. For this source, the same 
temperature profile as for IK Tau was used. 

The line profiles produced by our NH3 model calculations 
are shown in Figs. [Tl and I A. 1 1 overlaid on the measured spectra. 
The resulting NH3 abundances for the four observed stars range 
from 2 X 10"^ to 3x10"*' (see Table[T]i- While they are in Une with 
circumstellar NH3 abundances derived from the inversion lines 
alone and from IR spectroscopy, we note, as a caveat, that our 
model calculations did not consider the possibility of IR pump- 
ing of the lo - Oo transition. IR pumping of the H2O 1 10 - lo.i 
line in IRC-h10 216's CSE via various vi brational bands has been 
investigated by lAgundez & Cernichard (|2006), who found it to 
be the dominant source of excitation of this high critical den- 
sity line over much of the star's outer envelope. Consequently, 
their modeling suggests an order of magnitude lower H2O abun- 
dance than in voked ea rlier from SWAS and Odin observations 
dMelnick et al...200L.Hasegawa et al...2006.) . 



4. Summary and outlook 

The high critical density of the Jk - lo - Oo ortho-NH3 line 
allows investigations of NH3 in a new density regime of circum- 
stellar envelopes. However, the remarkably high abundances we 
determine for all our objects confirm and significantly strengthen 
the finding that this molecule exists in a variety of CSEs at levels 
not explained by current chemical models. Better constraints on 
the emitting regions will come from HIFI observations of a num- 
ber of NH3 lines, which are actually evident i n the (heavily spec - 
trally diluted) PACS spectrum of VY CMa dRover etal.ll2010h . 
Given their intensities, such observations appear eminently fea- 
sible. They would place tight constraints on future radiative 
transfer modeling that includes infrared excitation. Additionally, 
with its extremely wideband new generation digital correlator, 
the Expanded Very Large Array (EVLA) will allow simultane- 
ous imaging of many NH3 inversion lines. Nevertheless, given 



the weakness of these lines toward ordinary AGB stars, PPNe 
being an exception, EVLA observations will be challenging. 
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Appendix A: NH3 rotational and inversion line 
spectra and fits for IK Tau and IRC+10420 



Fig.|A3]shows, for IK Tau and IRC+10420, the HIFI spectra of 
the Jk = !() - Oo ortho-NH3 line and the spectra of the (J, K) - 
(1,1) and (2,2) para-NHa lines produced from our VLA data 
together with our best fit model predictions. 




50 50 

LSR Velocity (km s ) 

Fig.A.l. Spectra of NH3 transitions for IK Tau (left column) 
and IRC+10420 (right column). Top row: HIFI spectra of the 
Jk = lo - Oo ortho-NH3 line. Middle row: VLA spectra of the 
(J,K) = (1, 1) para-NH3 Une. Bottom row: VLA spectra of the 
(/, K) = (2, 2) para-NH3 line. The intensity scales apply for both 
sources. The lefthand ordinate gives main-beam brightness tem- 
perature in a 40" FWHM beam, while the righthand ordinate 
gives flux density in Jy units. In all panels the red dashed line 
represents our best-fit model prediction. The blue bars give the 
spacings and theoretical relative intensities of the main groups 
of hfs components dete rmined by Cazzo li et al. (2009) for the 
lo - Oo transition and bv lKukolichI (Il967h for the (1,1) and (2,2) 
inversion lines, respectively. 
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